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О трансформациях гравиполя с помощью задачи Алексидзе 

Результат сьемок гравиполя - значения 2(5) модуля градиента потенциала силь тяжести (МГПС- 
Т) на некоторой гладкой замкнутой поверхности 5, целиком охватьввающей все тяготеющие массьт 
ГП). Направлениє 2(5) неизвестно, позтому задачу трансформаций поля невозможно свести к линей- 
ньтм граничньм задачам. Приближенное однозначноє решение задачи можно получить на незамкну- 
той большой (17х1)? поверхности 5 . 

Один из путей перерасчета значений силь тяжести в региональном масштабе - решение нелиней- 
ной внешней граничной задачи Алексидзе для потенциала силь тяжести (в частном случає - для ура- 
внения Лапласа) 





Задачу (1) лингаризируют путем предположений относительно направления градиента потенциала на 
поверхности 5 (2). Зти предположения вносят погрешности, т.к. определение направления силь тя- 
жести едва ли не сложнее определения модуля силь тяжести. 

Проблема разрешимости и единственности задачи (1) принципиальна в гравиметрии: отвечаєт, 
можно ли точно однозначно определить силу тяжести во внешнем пространстве, зная точно силу тя- 
жести на замкнутой поверхности 5, целиком охватьвающей массьг и какая нужна дополнительная ин- 
формация. При одновременном определений значения и направления силь тяжести на поверхности 5 
вьіЧчИСЛениє силь  тяжести во внешнем пространстве есть решением задачи Дирихле. Но совмест- 
нье определения модуля и направления силь6 тяжести дороги и в массовой гравиразведке не прово- 
дятся. 

Решают задачу (1) разложением по неортогональньм функциям с помощью итераций 
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где и" - Кк -е приближение решения задачи (1) в плоском случає. 

В плоском случає определения потенциала, как гармоничной функции, если МГПСТ задан на за- 
мкнутом контуре, задача (1) не имеет единого решения. Если действительная часть и(х, у) голоморф- 
ной на сфере С, полностью охватьвающей 5, функции Й(2)с-и (х,у)Чіу (х,у) - решение граничной 
задачи, мнимая часть будет решением той же задачи. 

Показано, что задача (1) имеет единственноє решениє, если 5 есть многообразиєм Ляпунова и 

ке 
направлениє а; (и ) мало отличаєтся от некоторого заданного направления в смьгсле метрики І, (5). 

Наша цель - исследовать численную сходимость итераций и погрешность заменьт уравнения си- 
ль тяжести (ему удовлетворяют граничнье данньт6е задачи (1) на сфере С с радиусом Земли) уравне- 
нием Лапласа. Зта погрешность (разность решений соответствующих задач типа (1)) при размере 
площади грависьемок 100 км, варьирует в пределах 75 мГал «єх 51 мГал |3| при перерасчете значе- 
ний силь тяжести в свободном воздухе по формуле Пуассона. Максимальная погрешность перерас- 


чета во внешнем пространстве относительно области С равна тах -дри 2-а(а)/4, 2(а) - значениє 


силь тяжести на поверхности идеально сферической Земли. Предположение гармоничности силь 
тяжести в перерасчете на вьсоту радиуса Земли генерирует погрешность 245 Гал. 


Из-за расхождения векторов составляющих аномалий силь тяжести Ідя| ш |Я- |З 
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гут значительно отличаться от ТОЧНЬІХх. Если даже аномалия силь; тяжести 7-8 точна, ее нельзя 


пересчитать во внешнем пространстве: дифференциальное уравнение силь тяжести существенно за- 
висит от скорости изменения углов между направлениєем аномалии силь тяжести и координатньми 
осями. В локальньх аномалиях зто направление может варьировать, и всякоеє предположение о его 
поведений существено искажаєт перерасчет аномалий силь тяжести. Позтому важно, что Я -|я| со- 


впадаєт с большой точностью с проекцией вектора аномалий сильг тяжести на направление вектора 


де 
Численное решение задачи (1) для 2- и 3-мерного случая для проверки сходимости итераций (2) 
на модельньх задачах из (1| с известньми точньми решениями (ради апробации программ) в горизо- 
нтальной плоскости 5 на равномерной сетке с шагом 10 км доказало разрешимость задачи (1). 
Первое приближение направляющих косинусов є, - грубье приближения 0, 0, І, в ином вариан- 


те - 0.8, 0, 0.6. Число итераций не превьшало 10. Сравнениєе точньх и приближенньх значений МГ - 


ПСТ обнаруживаєт, что их граничнье значения во всех точках совпадают с точностью до 107", изме- 
няясь в пределах 0,07- 0,133, но они сами совпадают с куда меньшей точностью из-за плохой обусло- 
вленности задач с производньми в граничньх условиях (учитьваєм зто при решений практических 
задач). 

Решение с помощью итераций (2) внутренней граничной задачи 
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О) 


с началом координат (0,0,0) на расстояниим 0,2 от нижней кромки единичного куба С (его внешняя 
нормаль имеет 6 разньх направлений, 24 вспомогательньєе точки взять на кубе стороной 14, а а; 


равно 0, 0, 1) виявило следующую особенность. После 5-й итерации граничньєе значения удовлетво- 
реньг (погрешность Х196), но решение и его производньве не имеют ничего общего с точньм решени- 


10 за 
ем Пс и се производньчаи. Увеличениє числа итерации до 24 не увеличило точность 


а хі за 
приближения. Прибавка 25-го фундаментального решения, совпадающего с точньм, решило пробле- 
му: после 4-й итерации и граничньге условия, и значение градиентов и их производньх удовлетворе- 
ньг с точностью 107. 





Сходимость итераций (2) изучена и 
путем решения внешней задачи (1) для 
единичного круга с центром в начале 

5 координат (точное решениє 


ди ні 2 2 
з - дж изіпаху бХх), 





сі з сов(г - ф), с 8іп(т-нф), г - центра- 
ЛЬНЬЙ угол, Ф? - мера возмущения на- 


чального направления искомого гради- 
ента) с обрьвом итераций по условию 


ж - етади" | - І «107. Для сходимости 





дх, із | итераций (2) необходимо не только, 
вгаан чтобью точноєе решение удовлетворяло 
условию |1| 





іо (М,1)» 0, (4) 
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но и вьтбор такого ад, для которого в (2) не возникали бь козффициенть й, , не обеспечивающиєе 
вьшолнения условия (4). Условие (4) значит, что направление внешней нормали М неортогонально 
направлению І косой производной, т.е. | ни в одной точке границь 5 не совпадаєт с касательньм на- 
правлением. Условие вьшолняєтся для любой конечной части плоскости 5 для поля точечного источ- 
ника, а на всей бесконечной плоскости оно невьшолнимо. 

При расчете вместо потенциала У/ силь  тяжести потенциала Ї/ аномалиий силь тяжести ( 
У «СУ, где У - известньй потенциал нормальной силь тяжести) в линейньх приближениях гра- 
ничньх задач грави- и магнитометриий необходимо в качестве граничной функции у (5) взять анома- 


лию силь6 тяжести. В нелинейньх задачах переход к аномалиий изменяєт левую часть граничньх ус- 
ловий (1), ибо операция вьчисления МГПСТ не дистрибутивна: (гад МИ| ее вгад 2 --|огад М. Если 


у (5) - МГПСТ У/, определяем потенциал аномалий Ї/ из решения нелинейной граничной задачи 
Ли(х, ) з 0, х, є б, ці ЗМ, СІНО НО? ди, ЗУ, УЮ, )- у"(5)- |тадУЇ., (5) 


граничнье значения которой при нулевом потенциале нормальной силь тяжести совпадают с гранич- 
ньтми условиями (1). Вняснение численной скорости сходимости задачи (5) имеет отдельнь6й практи- 
ческий интерес. 
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